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Abstract 

In this paper we study the asymptotic behaviour of a nonlocal non- 
linear parabolic equation governed by a parameter. After giving the exis- 
tence of unique branch of solutions composed by stable solutions in sta- 
tionary case, we gives for the parabolic problem L°° estimates of solution 
based on using the Moser iterations and existence of global attractor. We 
finish our study by the issue of asymptotic behaviour in some cases when 
t oo. 



1 Introduction 

The non-local issues are important in studying the behavior of certain physi- 
cal phenomena and population dynamics. A major difficulty in studying these 
problems often lie in the absence of well-known properties as maximum princi- 
ple, regularity and properties of Lyapunov (see [5], [5]) and also the difficulty 
to characterize and determine the stationary solutions associated thus making 
study the asymptotic behavior of these solutions very difficult. 
In this paper we study the solution u(t, x) to the nonlocal equation 



(1) 



In the above problem uq and / are such that 

u eL 2 (n), /eL 2 (0,T,L 2 (fi)), (2) 
with T a arbitrary positive number, a is a continuous function such that 

3m, M such that < m < a(e) < M Ve G R. (3) 
The nonlocal functional l r is defined such that 

I P (.)(a;) : L 2 (n) -> K, u -> l r (u(t))(x) = J g(y)u(t, y)dy. (4) 

nns(x,r) 

Here B(x,r) is the closed ball of R" with radius r and g G L 2 (f2).. It is some- 
times possible to consider g more generally, especially when one is interested in 
the study of stationary solutions of (see [3]). 
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In physical point of view problem (JT|) gives many applications especially 
where g = 1 in population dynamics. Indeed, in this situation u may represent 
a population density and l r (u) the total mass of the subdomain Q n B{x, r) of 
f2. Hence (TTJ) can describe the evolution of a population whose diffusion velocity 
depends on the total mass of a subdomain of ft. For more details of modeli- 
sation we refer the reader to [7|. This type of equation can be applied more 
generally to other models including the study of propagation of mutant gene 
(see [H] , [H] , Q3] ) ■ A very recent study of this propagation was made by Ben- 
dahmane and Sepulveda [¥] in which they analyze using a finite volume scheme 
adapted, the transmission of this gene through 3 types of people: susceptible, 
infected and recovered. 

In mathematical point of view, when r — d where d is the diameter of f2 
problem (|T]| has been studied in various forms (see [S] . [5] . [5] . [T5] ) . 

However when < r < d, several questions from the theory of bifurca- 
tions have arisen concerning the structure of stationary solutions including the 
existence of a principle of comparison of different solutions depending on the 
parameter r and the existence of branches (local and global) of solutions. A 
large majority of these issues has been resolved in [3J. It shows that when a is 
decreasing the existence of a unique global branch of solutions and existence of 
branch of solutions that are purely local. Some questions may then arise: 

(i) The unique branch described in [3] it is composed of stable solutions? 

(ii) What about stability properties of the corresponding parabolic problem? 

The plan for this work is the following. In section 2 we give some existence and 
uniqueness results. Section 3 is devoted to stationary problem corresponding to 
(H|). In particular, we study in a radial generalisation of Chipot-Lovat 

results about determination of the number of solutions. We also establish that 
the unique global branch of solutions described in [3J is composed by stable 
solutions ftheorem l3.8[) . In section 4 firstly we address a L°° estimate taking to 
account LP estimate based on Moser iterations. Secondly we prove existence of 
absorbing set in LL\ , which allows us to prove the existence of a global attractor 
associated to ([T]) (see remark[5]) . Finally we obtain a result of stability properties 
of the corresponding parabolic problem. 

2 Existence and uniqueness results 

In this section we show a result of existence . We set V = Hq(H) and V its 
dual, we take the norm in V, \\-\\v such that 



< . , . > means the duality bracket of V and V. 
Then we have 

Theorem 2.1. Let T > 0, / € L 2 (0,T, V) and u e L 2 (£l), we assume that a 
is a continuous function and the assumption |3J) checked then for every r fixed, 
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r £ [0, diam(Q)], there exists a function u such that 
ueL 2 (0,T,V), u t e L 2 (0,T,V) 



it(0, .) — uq in f2 



f t (u,(f>) + J n a(l r (u(t)))VuV<t>dx =< f,cf>> in D'(0,T) V</> e H^il). 



(5) 

Moreover if a is locally Lipschitz i. e 

Vc 37 c such that \a(e) — a(e')| < j c \e — e'\ Ve, e' € [— c, c], (6) 
t/ien </ie solution of |3P is unique. 

Remark 1. Before to do the proof it is necessary to see that for r = problem 
|2P zs linear and the proof follows a well-known result (see J10\/). it is even 
when r — diam(fl) (see flf). We will focus therefore in the following where 
r s]0, diam(Q)[. 

Proof. For the existence proof we will use the Schauder fixed point theorem. 
Let w £ L 2 (0,T,L 2 (n)) we get 

t—>lr(w(t)), 

is measurable as a is continuous then 

t — >■ a(l r (w(t))), 
is too. The problem of finding u = u(t, x) solution of 

u e L 2 {0, T, V) n C([0, T], L 2 (fi)) ut £ L 2 (0, T, V) 
u(0, .) = u 

jt(u,4>) + J n a(l r (w(t)))VuV(l>dx=<f,(f>> in D'(0,T) e H%{Sl), 

(7) 

is linear, besides ([7]) admits a unique solution it = F r (w) (see[TU],[7]). Thus we 
show that the application 

w — >f r H=u, (8) 

admits a fixed point. Taking w = u in (J7J| we get using and the Cauchy- 
Schwarz inequality 

~l«l2 + "»||u||^<|/UHv, (9) 
H-llv is the usual norm in V and |/|* is the dual norm of /. We take 

T i 
|2 



Ml2(o,T,V) = \ / IMIv* 







Using Young's inequality to the right-hand side of it follows that 
By integrating (fTQ)) on (0, t) for i < Twe obtain 



t 



\\u{t)\l + | J \\u\\ 2 v dt < ±\u \ 2 2 + ±.J |/| 2 . (11) 
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We deduce that there exists a constant C = C(m, uq, f) such that 

Hmo,T,v) < C (12) 

Moreover 

< u t , v > + < -div(a(l r (u(t)))Vu), v >=< f,v> Vv € V, 
This gives us 

|ut|*<M||u||vr+|/U. (13) 
By raising (|13p squared and using the Young inequality we have that 

\u t \l<2M 2 \\u\\ 2 v + 2\f\l (14) 

By integrating (fT4"]l on (0,t) and assuming (IT21 we obtain 

Ml 2 (o,t,v") < C, (15) 
with C" = C'(m, M, /, uo) and C" is independent to iu. It follows from (fT2|) and 

(113 

l w tL 2 (o,T,y) + Ml 2 (o,t,v) ^ -K) (16) 
with R = C 2 + C' 2 . From (fT2l and the Poincare inequality it follows that 

|w|i,2(o,r,L 2 (n)) < R\ (17) 

By setting 

R 1 =max(R',R), (18) 

and associating (fTT)) and (fT8"]) . it follows that the application F maps the ball 
5(0, of L 2 (0, T, L 2 {Vt)) into itself. Moreover the balls of H 1 ^, T, V, V) are 
relatively compact in L 2 (0, T, L 2 (£i)) (see (TUj for more details), (fTB|) clearly 
shows us that F(B(0, Ri) is relatively compact in 5(0, i?i) with 

5(0, JZi) = e L 2 (0, T, L 2 (r!)); |«| La(0 ,T,L»(n)) < 

In order to apply the Schauder fixed point theorem, as announced, we just need 
to show that F is continuous from 5(0, i?i) to itself. This is actually the case 
and completes the proof of existence. 

We will now discuss the uniqueness assuming of course that assumption (jB]) 
be verified. Consider u\ and u 2 two solutions ([5]), by subtracting one obtains in 
D'(Q,T) 

d f 

— {u x -u 2 ,v)+ I (a(l r (u 1 (t))Vu 1 (t)-a{l r (u 2 (t)))\7u 2 (t))\7(j)dx = £ H^(Q). 

Q 

(19) 

Since 

o(/ r (ui(t)))Vui - a{l r {u 2 (t)))Vu 2 {t) = (a(l r (ui(t))) - a(l r (u 2 (t)))Vui(t) 

+ o(l P («2(t)))V(ui(t)-«a(*)), 

(20) 
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we get 

d f 

— (m - u 2 ,v)+ / a(l r (u 2 (t)))V (ui(t) - u 2 (t))V (j>dx 
n 

= - J {a{l r {ui(t))) - a(l r (u 2 {t)))\7 Ul \7<j)dx V<j> G H^(Q). 
n 

(21) 

Moreover ui,u 2 G C([0, T], L 2 (Q)) there exist z > such that 

l r (ui(t)),l r (u 2 (t)) e[-z,z}. (22) 
Taking v = u\ — u 2 in (I2ip . it comes easily by Cauchy-Schwartz inequality and 

-U2I2 +m||tti -tta||v < 7Rr(«i(*)) - Jr("2(*))|||«i||vil«i - w 2 II v- (23) 
We get in [3J 

\l r (u(t))<C\B(x,r)nn\^s\g\ 2 \u(t)\ 2 < |s| 2 Kt)| 2 , (24) 

where C a constant, |0| represents the measure of f2 and n V 3 the maximum 
between the dimension n of f2 and 3. By using (|24l) . (|23|) and the Young in- 
equality 

1 , 9 m , 
~ 2m 2 

We deduce 

— |ui - U2I2 + "i||m - u 2 ||y < p(t)|«i - u 2 | 2 : (25) 

with 

p {t) = I( 7 c|0|^3| 9 | 2 || Ul || y ) 2 e ^(CT), 



which leads to 



-H2I2 < p(t)\ui -u 2 \l- (26) 



Multiplying ([H) by e" /<? p( s ) ds it follows that 

e -/>( s )rf s ^| Ui _ U2 |2 _ p ( t ) e -/oP( s ) ds | Ml - U2 |2 < 0. (27) 

Hence 

|.{ e -/o*PW<i«| Ul _ M2 |2}< . (28) 

t 

— J p(s)ds 

This shows that t 1 — > e \u± — u 2 \ 2 is nonincreasing. Since for t = 0, 

ui(0, .) = « 2 (0, .) = u . 

This function vanishes at and nonnegative, we conclude that it is identically 
zero. This concludes the proof. □ 
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3 Stationary solutions 

Consider the weak formulation to the stationary problem associated to ([IJ 

(p ,(-div(a(l r (u))Vu) = f dans Q 

1 r, \ueH^(Q). [ ' 

3.1 The case r = d 

By taking (f> the weak solution of the problem 



-A(j> — f dans il 



(30) 



we get due to a Chipot-Lovat [8] results that 

Theorem 3.1. Let a be a mapping from K into (0, oo). TTie problem (Pd) has 
many solutions as the problem in 1R 

Ha(n) = l d (4>), (31) 

with fi = l ( i{u d ). 

Remark 2. Theorem \3.1\ allows us to see where a is increasing that the problem 
Pd admits a unique solution and determine for a given a the exact number of 
solutions (Pd)- However it is difficult or impossible to adapt the proof of the 
theore m] 3. 1\ in case < r < d. 

3.2 The case < r < d 

As announced in the introduction we focus our study to the case of radial 
solutions of {Pd)- We will assume £1 is the open ball of K n with radius d/2 
centered at zero. We set 

L 2 r (Q) = {u e L 2 (n) 3u e L 2 (]0, d/2[) such that u(x) = u(||a:||)}, 

and we also assume that 

/ G Ll(Q) 
g E L 2 r (Q) 

a g W 1 '°°(M.) 1 inf a > (32) 

/ > a.e in fi 
g > a.e in f2. 

We start by giving in some sense in a linear case a result that will be used 
later to explain the asymptotic behavior. 

Proposition 3.2. LetA,B £ C(£l) be positive radial functions such that A < B 
in fl and also f,h £ L 2 (fl) two positive radial functions. Let u G Hq(£1) the 
radial solution to 

- div(A{x)Vu) = f in ft, (33) 

and 

- div(B(x)Vu) =h in tl. (34) 

Then f < h a.e in fl. 



G 



Proof. We proved in |3J that if u is a the radial solution of (f3"3")l then for a.e t in 

[0,d/2], 

t 

u'(t) = ~J^J( S -T- 1 f( S )ds. (35) 

From ([55]). (|MD and we obtain 

t t 

iJ(t) 



□ 



i(/j/( -r/(^ s =/(-rM^ s . 



Since A < B in and /, h > with / ^ 0, /i ^ hence / < g. 



In a nonlocal case, some results of existence of radial solutions and compar- 
ison principle between u r , Ud and uq has been demonstrated in [3 . It is also 
proved that if we set for all r G [0, d] 

I r := [inf l r (</>), sup (36) 

Here </> denotes the solution of 

(-A<j) = f dans n 

By the inclusion or not of I r at an interval of R we somehow generalize the 
theorem l3.ll 

Lemma 3.3. Let r G [0, d]. Assume that hS'A) holds true and there exist < 
mi < rri2 such that 

a(mi) = max a 0(1112) = min a (38) 

[mi,m 2 ] [1111,1112] 

I r C [mifl(mi), m 2 a(m 2 )]. (39) 
TTien (P r ^ admits a radial solution u and 

mi < Z r (u) < m2 a.e in fi. (40) 

For the proof, we refer the reader to [3]. 

Generalizing this construction type of the diffusion coefficient a we obtain 

Proposition 3.4. Let r G [0,d\. Assume that \3 6 A) holds true and there exist 
an odd integer n\ and n\ + 1 positive real numbers {?7ij}i = o...n 1 , with too = 
and for all i G {0, . . . , n% — 1} we /iawe mj < m !+ i. Moreover 

a(rrii) = max a; a(m 1+ i) = min a Vi € {0, 2, . . . , n\ — 3, Jii — 1} 

[mi,mi+i] [mi,mi-(-i] 

J r C Pi [mia(mi),m l+1 a(m i+1 )} 

i— 0,2,...,m— 3,ni — 1 

(41) 

TTiew (P rv ) admits at least "'^ radial solutions {ui} ie ^ ,2,. ..,711—1} such that 
mi < l r (ui) < irii + i Vi G {0, 2, . . . , ni — 3, ri\ — 1}. 
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Proof. The proof here is by induction. Indeed we set 

Vm — { If condition BT]) is satisfied then(P r ) admits at least 1 ^~ — solutions.} 

By using lemma [3~l?l with mi = and m.2 = mi, it is easy to prove for n\ = 1 
that 7 ? „ 1 is true. For ni > 1, This procedure can be repeated to prove that if 
7 ? ni-2 holds true then V ni holds too. □ 



Example 1. Let us see a function a satisfying proposition \3.4\ For this, we 
consider the case n\ = 3 and r G (0, d] . Considering V3 6 2\) and the strong 
maximum principle we get minl r > 0. Taking 

maxl r a(0) 

mi:=2 w a(mi):= — 

with a(0) > and also a decreasing on [0, mi] then we prove lemma \3~3^ condi- 
tions. 

By repeating this process with m,2 > mi and setting 

minl r maxl r 
a(m 2 ) := , m 3 := 2 ■ 



m 2 a(m.2) 

with a(ffi3) := a ^" 2 ^ and also a is decreasing on [m2,ms\- This shows the 
existence of such a. 




m 2 m 3 

Figure 1: The case n\ = 3 

In the representation of a we have deliberately left, on solid line parts of 
the curve satisfying the conditions of proposition ^. 4\ and dotted line one without 
constraints. This situation are explain in the figure^ 
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Remark 3. As previously announced, the proposition \3.4\ generalizes a certain 
point of view Theoren A3.1\ However it does not accurately determine the exact 
number of solutions of (P r ) and the bifurcation points of branch of solutions. 
We have shown in way to solve this problem by using the linearized problem, 
the principle of comparisons obtained in and the Krein-Rutman theorem. 

3.3 Stable solutions of (P r ) 

Definition 3.5. Given a domain Q. C W 1 , a solution u r G H^fl) of (P r ) is 
stable if: 

G ffo(fi) G Ur (0) := / a{l r {u r ))\\7(j)\ 2 - / a' ' {l r {u r ))l r {<j))Vu r V '(/> > 0. 



(42) 

Definition 3.6. Given u : [0,d] — > Hq(CI), the graph of u is called a (global) 
branch of solutions if 

(i) ue C([0,d], H£ (SI)), 

(ii) u(r) is solution to (P r ) for all r in [0,d]. 

u is called a local branch if it's defined only on a subinterval of [0, d] with positive 
measure. 

Before concluding this section, we will focus into the case a nonincreasing to 
prove the stability of the global branch of solutions. Assume for all r G [0, d], 
u r is a solution to (P r ) and 

< l r (u r )(x) < fid for a.e x 6 ft. (43) 

Assume that there exists a solution fid to (|3"Tj) such that 

a(iJLd) = niin a and a(0) = max a. (44) 

We prove in [3] 

Theorem 3.7. Assume A Sty) , j43\ ), j44\ ) an d fHP holds. Assume in addition 
that a G W ,oc (M) and for some positive constant e, it holds that 



where C\ is a constant dependent to il. Then 
(i) For all r in [0,d], (P r ) possesses a unique radial solution u r in [uQ,Ud\; 
(ii) {(r, u r ) : r G [0,<i]} is a branch of solutions without bifurcation point; 
(Hi) it is only global branch of solutions; 

(iv) if in addition, a is nonincreasing on [0,fid] then r M> u r is nondecreasing. 
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Remark 4. It is very difficult to obtain property (iv) for any a. However when 
a is nonincreasing provide us important information for studying the stability 
of this branch of solutions. 

Corollary 3.8. Let u\ the smallest solution to (Pd)- Assume J3jj|) and 
holds true and there exists a solution fid to VS1\) satisfied {44\ l- Assume in 
addition that a £ W 1 ' 00 ^), u\ satisfied |^ff| ) and for some positive constant e, 
it holds that 

1 

where C\ is a constant dependent to CI. 

Then {(r,u r ) : r £ [0,d]} is the only global branch of solutions starting to 



Cilffbl/blo'looJ-^+^TTT^T < !. ( 46 ) 



Proof. The fact that {(r,u r ) : r £ [0, d]} is the only global branch of solutions 
results from theorem l3.7l We will now show that this unique branch of solutions 
is stable and start at r = d by u d . For this we consider without loss of generality 
(Pd) admits two solutions u d and u 2 d such that u l d < u% We denote by /ii and 
/i2 respectively solutions of (f5T|) corresponding to u d and u\ (see figure[3]). It is 
easy to see that /ii and fi2 satisfied (|4"4"|) . 

Assume {(r,u r ) : r £ [0,d]} is the only global branch of solutions starting 
to u% Then we get C 1 \g\ 2 \f\2\a'\ O0 ^- <i ^ 2+( ]- a j^i < 1. In this case, using 
theorem !3.7l we get (P r ) possesses a unique radial solution u r in [i*o,it^] and the 
mapping r h-> u r is nondecreasing. By continuity of this mapping, we can find a 
ro £]0,d[ such that u ro = u l d for a.e x £ CI. This means that u d is a solution 
of (P ro ) . This gives us an absurdity and concludes the proof. □ 

We are now able to prove: 

Proposition 3.9. Under assumptions and notation of corollar u\3.8\ the global 
branch of solutions described in theore m\3. 71 is composed by stable solutions. 

Proof. For all r £ [0,d\, let u r be a solution belonging to the global branch of 
solutions described in theorem !3.7l By using the linearized problem of (P r ), we 
get £ H%(Cl) 



a(l r (u r ))\V(j)\ 2 - / a'{l r {u r ))l r ((j))Vu r V(j) > 



n si 



(47) 



Taking into account that |Vu r |2 < ^(^) |^f"^7w^ T) wnere ^(Cl) designed the 



si 

Poincare Sobolev constant. We obtain 



inf o(/r(«r))|V^| 2 — ^ |ff|2|a'|oo,[- e , Ml +e] | Vli^ 1 2 | V0|| . 

I/I2 
a(l r (u r )) 

a(l r {u r ))\V(j)\ 2 - J a'(l r (u r ))l r (<j))\7u r V<f> > 
h 

|V0| 2 (inf a(l r (u r )) - d LgHalod-e^+e 



I/I2 



inf a(l r (u r )) 

SI 



(48) 
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Moreover by assumptions (j^51) and (pPtj) we get a(fj,\) < inf a(l r (u r )). 
Thus (|46|) becomes 

Cl|ff| 2 [/l2la / loo,[-e, At<i +e] , a < 1. 

mi a t r (w. r 

Q 

We deduces 

y a(/ r (u r ))| V0| 2 - y a' {l r {u r ))l r {4>)Vu r V (f) > 0. 
n n 

This concluded the proof. 



aim 

a(\mu_l) 



Figure 2: case of 2 solutions 





"i 



r=0 



r=d 



Figure 3: Branch of solutions 
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4 Parabolic problem 

4.1 L°° estimate 

In what follows we obtain L°° estimate of the solution ([T]) from L q estimate. 
The method we use is based on iterations Moser type, for more details on the 
method see |14| . 
We get 

Theorem 4.1. Let n > 3 and u a classical solution of (QP defined on [0, T). 
Assume that p > 1 and q > 1 smc/i i/iai ^ + | = 1. Suppose further that 
U q = sup t<T \u(t)\ q < oo, f e i°°(0, cx), L«(fi)). 7/p < ^ </len U °° < 00 • 

To prove this theorem we need the following proposition: 



Lemma 4.2. Consider u a classical solution of (QJ) on [0, T), r > 1 and p > 
1 sitc/i t/iat - + - = 
supt<T\u(t)\ r } and let 



1 smc/i that i + | = 1 p < ^rj. We take U r = max{l, |ito|ooj U r 



p{n + 2) 

cr(r) - 



2[r(2p — pn + n) + np] 
Then there exists a constant Ci — C'2(Q,m) such that 

U2r<[C 2 \\f\\L-(0,ooMm)r {r) r a{r) Ur. 

Proof. Multiplying (Q]) by u 21 "" 1 and then using the Holder inequality yields 

~ J u 2r dx + m^-j± J \V(u r )\ 2 dx < \f\ q \u 2r -\. (51) 
a n 

As 

\u 2r -\=\vr\^, (52) 
by taking w = u r in (f5"Tj) and ([5^)1 . we get easily 

^^Ha + m^i^|VHl<|/| 9 K p , (53) 



with a = . Let 5 such that 



£-0+^. (54) 

with 2* = We claim that /3 e (0, 1). 
In fact 

2nr - (n- 2)(2r - l)p 
^ ~ (n + 2)(2r- l)p ' 

Sincep < ^ then /? > 0. As well as (ra+2)(2r-l)p > 2w-(n-2)(2r-l)p 
implies that j3 < 1 this prove that /? G (0, 1). 

Using an interpolation inequality (see [14]) in (j53|) and (|54|) . we get 

1 d ' .12 , 2r-U a , ,1-/3 



2r -^H2 + m -^-l v H2< KHir p ) . (55) 
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Applying Sobolev injections in (|55[) . we have 



Id, ,n 2r — 1 ._ l9 



and also 

1^ii9 2r — 1 ._ l9 



l/l,(-Y ' \< a c^ a 

1 m y 



(56) 



2r /' 



(57) 

Since /3 G (0, 1) and f G (0, 1) it is clear that a(1 2 p) G (0, 1). Applying Young's 

Afe obtain 

a(l-/3) 



inequality to ([57]) with + 1 - = 1. We obtain 



l^ii? 2r — 1 ._ l9 



ir^ 



with (5=1 — 



"(l-ff) 



Joining the fact that a ^ 2 ^ G (0, 1) and <5 < 1 to (f58")l . we deduce 



1 d 



3r - 2 



+ ™^2-\vw\i < |/|| ( - ) H7" <?* 



2r 



»(l-<3) 



We set 



2r 2 



"(1-/3) Pa 
2m(r) - 1 = 1 ^ and 2p(r) = 



m \ I,-, 

27j |v ^ 

(58) 



(59) 



(JSHJ) becomes 
1 d 



3r - 2 



27d7 H ^ +W: 2r 2 



|VHI<I/II(- 

m 



2ra(r)-l 



\w\\ p{r 'CK (60) 



This gives us taking into account that r > 1 



2r\ 



2ra r (r) 



j t \w\i + m\Vw\i < 1^1 Mi pW mC* . 
By a calculation we can verify that 

2nr -(n- 2)(2r - l)p 



(61) 



p(r) 



2r(p(n + 2) + n) - 2n(2r - l)p' 



and also that p(r) G (0, 1). 

Using the Poincare Sobolev inequality and that p(r) < 1 in (f6"Tj) . yields 



I /2r\ 



where Ci (SI) designed the Poincare Sobolev constant. Noticing that 



(62) 



dt 



Ci(Jl) 



2r\ 



2ru(r) 



w\i<\f\l[-J KmC. 

(63) 
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and integrating ([55)1 on [0,t) we get 

(~ \ 2r<r(r) 



(64) 



Since 

|to(0)|i = / «;(0) 2 dx = / u(0) 2r dx < \n\\u(o)\ 2 z < m? r , UK l 



(|6"4")l and ([B3|) gives us 

2r(j(r) 



XJll < \Q\U? + \\f\\l x{0tO2!LHn)) m C^U?. (66) 

Whereas i > 1, 2ra(r) > and cr(r) = ^3 it follows that 

with C2 = C2(il,m). This completes the proof of Lemma. □ 
We have also 

Lemma 4.3. Let r > 1, n > 3, p < ^ a7ld CT ( r ) = 2lr-(2 P - ( pn + +»)+np| </ieri we 
get 

a{2 k r) < 9 k a{r) Vfc G N, 

with 9 £ (0, 1). 

Proof. By asking ci = EilLtHl^ C2 — (2p — + n) and C3 = yields a(r) — 
rc ^ c with ci, C2, C3 G .By taking 9 = 1 — 2C2+C ^ ne P ro °f °f this lemma is 
deduced by reasoning by induction. □ 

Returning now to the proof of the theorem, by lemma B31 we get 

Proof. 

U 2r < [C 2 ||/IU»(o,oc,L«(n))]* (r V ff(r) tf r . 
By iterating this equation by taking r — h,r — 2h, r = 2 2 h, etc, we obtain 

U-zh+lfr — [Ci ||/||l°°(0,oo,L9 m r2^h^ Uh , 

with 

Ai := a(h) + a(2h) + o{2 2 h) + .. + a(2 k - 1 h) + <r{2 k r) 

et 

A 2 := a(2h) + 2a(2 2 h) + 3a(2 3 h) + ... + (k - l)a{2 k - 1 h) + ka{2 k r). 

To complete the proof we just need to show that Ai,A2 < +00. Indeed by 
leinma POl 

k 00 / , n 

Ai < V a^aih) < V aPath) = , \ < 00. 
z — ' ^— ' (1 — a) 

Noting also that 

a(2 k h) < 9 k - 1 a{2h) Vfc G N*, 
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it follows 

A 2 < m^v&h) < jr fia'"- 1 a(2h) = < oo. 

fj,= l n=i ^ ' 

This completes the proof of the theorem. □ 



4.2 Uniform estimate in time 

We prove in what follows an estimate of u in i°°(M + , i/g (17)). We get 

Theorem 4.4. Assume that f e £ 2 (17), g £ H 1 ^), uq £ Hq(Q) and a 6 
W 1,oc (R) withmi a > 0. Then a solution u of is such thatuE L°°(R+, H%(Q)). 

]R 

Taking a spectral basis related to the Laplace operator in the Galerkin ap- 
proximation (see |16|) we find that — Au can be regarded as test function in 
L 2 (0,T,L 2 (17)) for all T > 0. By multiplying (D by -Au(t) and integrating 
over 17 

(u t , -Au) + {-div(a(l r (u))Vu), -Au) = (/, -Am), (68) 

and also 

^Nlv + (-a(i r (iO)Au, -Au) + (-a'(l r (u))Vl r (u).Vu, -Au) = (/, -Au). 

(69) 

Here (., .) is the usual scalar product on L 2 (17). Taking to account 

\Vl r (u)\ 2 <K\\g\\ Hl{n) \Vu\ 2 , (70) 
where if is a constant depending of 17. It comes 

|(-o / (/ r (tt))Vi r («).V«,-Au)| < iflMlffi^la'UMIvlAula (71) 
Now from (|7ip and ([6"9"| we have 

i||| U || 2 y + m|A U | 2 -X||.g|| ffl(0) | fl '| 0O |l«lly|Au| 2 < |/| 2 |A U | 2 . (72) 

By using Young's inequality ab < ^a 2 + ^b 2 , we get 

1 . .,o 1 

,.M\v< 

at mm 
In order to apply the uniform Gronwall lemma to (J73J) we start with a small 
estimate. Remember that 



\u\\ 2 v < -|/| 2 + -{K \\g\\& (a) ?\\a%U\ut- (73) 



^H2+"*||«||v<^l/li (74) 

where A is the principal eigenvalue of the Laplacian operator with Dirichlet 
boundary conditions. 

By integrating on [t,to) we get 

t+t t+t 

\u(t + t Q )\ 2 +m f \\u\\ 2 v ds< f A-\f\ 2 ds+\u(t)\l (75) 



m 
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and also 

t+t 



[ 

J A to 2 ' J lz to 

t 



u\\ 2 v ds<-^\f\ 2 ds + -\u(t)\ 2 . (76) 



Let po > such that |u(t)|| < p\- By setting 

ox = l Cl (0)VlL«3 aa = ^|/|i a3 = ^|/|l+-pg, 

m TO TO TO 

we obtain by using uniform Gronwall lemma to (|73[) 

Ki + *o)||v < + a 2 )exp(a 1 ) Vt > 0, i > 0. (77) 
to 

Hence u £ L°°(to, +oo, J/q (17)). By using ([T5|) and the classical Gronwall lemma 
it is easy to see that u £ L°°(0, fo> -Ho (^))- This completes the proof of the 
theorem. 

Remark 5. This theorem show us the existence of absorbing set in Hq(Q,). By 
considering S(t) the semigroup associated to the equation (QJ) defined by 

S(t) : L 2 (17) -> L 2 (n) 
u Q (->■ u(t), 

wii/i u(t) a solution of ([!]). .As a result of the theore m\4-4\ an d the compact 
embedding of Hq(CI) into L 2 (17) we deduce that the semigroup S(t) possesses 
a global attractor. Indeed it is easy to show the existence of absorbing set in 
L 2 (Q), the main difficulty here is to show that S(t) is such that 

VB C L 2 (17) bounded, 3t = t (B) 

such that |)S'(t)B is relatively compact in i 2 (17). C^) 

t>t 

This property known as S(t) is uniformly compact for t large can be proved by 
using theorem \4-4\ an d the compact embedding of Hq(£1) into L 2 (17). 



4.3 Asymptotic behaviour 

In this part we are interested in asymptotic behaviour of a weak solution of (fTJ) . 
Our main interest here is the radial solutions. By radial solutions we means 
u(t, \x\) = u(t, x). As in the stationary case 17 is a open ball of M™. Remember 
that 

L 2 (17) = {v e L 2 (17) 3v £ L 2 (}0, d/2[) such that v(x) = v(\\x\\)}. 

In order to not make confusion between uq the solution to (Po) an d the initial 
value of (jTJ) , we will take u° the initial value of |T|). 

Theorem 4.5. Assume that f,g € L 2 (17), a is a continuous function and the 
assumption |3J) checked then (J\l admits a radial solution. 

Proof. Let w G L 2 (0, t, L 2 (17)) it is clear that l r (w) is radial and also a(l r (w)). 
Thus by © F r maps L 2 (0, t, L 2 (17)) into itself. The proof now follows by using 
arguments similar to those used in thcorcm l2.ll □ 
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Assume now 

f,g>0 in n (79) 

and 

u <u° < u d , (80) 

with u° the initial value to |1| and uq and u d respectively the solution of (Po) 
and (P d ). 

We can now give a stability result assuming that ([T]) admits a unique solution. 

Theorem 4.6. Assume \ 19i and f,g G L^(fl). Let u, u d and uq respectively 
the solution of (Pd) and (Pd)- If 

uo < u° < u d , 

then 

uq < u < Ud Vf . 

Proof. Let 

S = {t\ l(u(s)) e[0,l d (u d )] Vs<t}. (81) 
It is easy to prove that S contains (see l80|) . By setting 

t* = sup{t \t G S}. (82) 
By continuity of the mapping t — > l d {u{t)), we get 

ld(u(t*)) G [0, ld{ud)]- (83) 

By using ^ and (P d ) we get in V(0,t*) 

j t {u d -u, (f>)+ J a(ld(u))V(u d -u)V(t> = - J (a(ld(u d ))-a(ld(u)))\7udS7<j) V</> G H%(Q). 
n n 

(84) 

Choising = (ud — u) , ([54]) becomes 

7^|(Mrf-ti)~l2+ J a(ld(ud))\S7(u d -uy\ 2 = J (a(l d (u))-a(l d (u d )))Vu d V(u d -u)~ . 
n n 

, . (85) 

Since a is nonincreasing (a(l d (u)) — a(l d (u d )) > Vi < t*) hence proposition l3.2l 
yields 

J {a(l d (u)) - a(l d (u d )))Vu d V(u d - u)~ < 0. (86) 
h 

Thus 

^|(wd-u) - |2 + a(^(wd))|V(u d -u)-|| <0 (87) 
Applying Poincarre Sobolev inequality we get 

\j t \{u d - u)-\l + C 2 \(u d - u)-\ 2 2 < 0, (88) 

this prove 

l{e 2t ^\(u d -u)-\l}<0. (89) 
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Moreover (u d — u) (0) = (ud — u°) = it follows that Ud> u Vt G [0, t*]. In 
the same way we can also prove uq < u Vt G [0, t*]. It follows 

mo < u < u d Vt e [0, t*] (90) 

To finish we just need to prove that t* is very large, this is typically the case. 
Indeed if t* < oo then 

l(u(t*))=0 or l d (u d ). (91) 
From (J75JI and ([50]) we deduce 

u(t*) = uq or u(t*) = u d . (92) 

Due to the uniqueness of we deduce that t = oo. This shows that 

Uq < U < Ud Vt, 

and achieve the proof. □ 

Remark 6. The fact that \u(i)\% is not a Lyapunov function that is to say de- 
creases in time, makes very complex the study of certain asymptotic properties 
of our problem. Indeed under our study it is tempting to show that for r fixed 

r e]o,d[ 

u(t)->ul in L 2 {n), 

where u is the solution of 0) and u], the solution belonging to the stable global 
branch described previously. A numerical study would be a great contribution to 
try to carry out some of our theoretical intuitions. 
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